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We study the behavior of an open quantum system, with an A''-dimensionaI space of states, 
whose density matrix evolves according to a non-unitary map defined in two steps: A unitary step, 
where the system evolves with an evolution operator obtained by quantizing a classically chaotic 
map (baker's and Harper's map are the two examples we consider). A non-unitary step where the 
evolution operator for the density matrix mimics the effect of diffusion in the semiclassical (large A'^) 
limit. The process of decoherence and the transition from quantum to classical behavior are analyzed 
. in detail by means of numerical and analitic tools. The existence of a regime where the entropy 

' grows with a rate which is independent of the strength of the diffusion coefficient is demonstrated. 

, The nature of the processes that determine the production of entropy is analyzed. 
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I. INTRODUCTION 

> 

CO I Decoherence has been recognized in recent years as one of the main ingredients needed to understand the origin of 
CO ■ the classical world from the fundamental quantum laws Decoherence is a process whose origin is conceptually 

simple: It is due to the entanglement between the system and its environment that is created in the course of their 
interaction. As a consequence, the environment keeps a record of the state of the system, that loses its purity. 
Only a small subset of all possible states of the system (the so-called pointer states) are relatively stable against the 
interaction. They are the ones which are less likely to become entangled with the environment. In the vast majority of 
cases, when the state is a superposition of pointer states, the information initially stored in the state of the system can 
never be restored since it irreversibly flows into the correlations with the environment. A basic question one should 
ask in this context is how fast does the information flow away from the system. This can be studied, for example, by 
' analyzing the evolution of the entropy obtained from the reduced density matrix of the system. This has been done for 
^ , a variety of cases (see |^ for a review) and it has been recognized that this process has unique features if the system 
has a classically chaotic counterpart. In fact, as conjectured in |^,^, the rate of entropy production in such cases has 
a regime that is independent of the strength of the coupling between the system and the environment and is entirely 
. ' determined by the dynamical parameters characterizing the chaotic evolution. This conjecture was analyzed in the 
^ , literature mostly using numerical tools P-pl| . In this paper we will present a study of this problem for some systems 
■ that are simple enough to enable both a rigorous numerical treatment and some analytic estimates. We consider here 
a quantum system with a finite dimensional Hilbert space {N is the number of dimensions, and we are interested in 
learning about the behavior of the system in the large N limit) . For such system we define an evolution operator for 
the density matrix in two steps as follows: first we consider a purely unitary evolution defining an operator U which 
is such that it corresponds (in the large N limit) to a classically chaotic map. Then, we define a non-unitary map 
for the density matrix in such a way that (again, in the large N limit) it mimics the effect of the interaction with an 
environment producing the same effects one observes for a Brownian particle (namely, diffusion). For such system, we 
developed a numerical code enabling us to efficiently evolve the density matrix and study in particular the evolution 
of the entropy. Our aim is to present solid numerical evidence supporting the conjecture presented in and, by 
combining the numerical work with analytic calculations, develop new intuition on the main processes contributing 
to entropy production for this kind of systems. 

The paper is organized as follows: In Section II we review some basic elements of the theory of quantum maps. 
We focus, as in the rest of the paper, on two specific examples: baker's map (the paradigmatic example of a fully 
chaotic system) and Harper's map (an example of the wide class of kicked maps with mixed phase space). There are 
no new results presented in this section and the reader with experience in the theory of quantum maps can easily 
skip it. In Section III we describe in detail the model for decoherence that we study in this paper. In Section IV we 
present the main results concerning the behavior of the entropy as a function of time and the evolution of the Wigner 
function. Finally, in Section IV we present our conclusions. An Appendix contains technical details about the phase 
space representation we use in this paper (the discrete Wigner function). 
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II. QUANTUM MAPS 



The construction of the quantum analogue of a classical map follows two well defined steps : a kinematical one, 
where the nature of the Hilbert space is defined in relationship with a specific phase space structure, and a dynamical 
one, where a unitary operator defines the evolution for a finite time step. For the study of chaotic behaviour, a finite 
phase space is required and therefore the question of boundary conditions arise. The simplest case is the torus, where 
periodic boundary conditions are assumed for both the coordinate and momentum representations. The most general 
quasi-periodic boundary conditions are 

(g + l|7A)=e»2-'^'(g|^) (1) 
(p+l|V-) = e-'2-x^(p|^) (2) 

where Xq Xp ^-re fixed, arbitrary real numbers between and 1 (27rxg and 27rxp are called Floquet angles). These 
conditions result in a finite dimensional Hilbert space. This space's dimension, N , is related to Ti by the relation 

2TihN = 1 (3) 

which signifies that phase space (of area equal to unity ) is spanned by N states (of area 2T:h). The position and 
momemntum eigenvalues in this finite dimensional space are 

kn) = |^^i^), n-0,l,...,7V-l, (4) 
\Prn)^\^^^), m = 0,l,...,iV-l. (5) 
These position and momentum eigenstates are related by a discrete Fourier transform: 

{pM = ^e-^(™+>^')("+>^^) ^ (G^^'^l™ (6) 



The values of and Xp specify different Hilbert spaces. In the present paper we will use Xg = Xp = ^, corresponding 
to antiperiodic boundary conditions. Cyclic shifts on these two bases are implemented by unitary operators 

l^\qn) = k„+i) (7) 

"l^bm) = \Pm+l)- (8) 

These operators and their powers can be combined to produce unitary displacement operators in phase space 

LI(Ap, Aq) = Z^Ag^Apg.^ApAg ^g) 

for integer values of Ap, Ap. These are the analogues of the Weyl displacement operators in the continuous case. 
Phase space in this context is then assimilated to a discrete N x N grid which will be useful for the representation of 
quantum effects in phase space (see Appendix). 



A. The baker's map 

The baker's map is one of the simplest systems displaying strongly chaotic behaviour. In spite of its simplicity, it 
has a very rich dynamical behaviour, both in its classical and its quantum versions. The map is an area preserving 
transformation defined in the [0, 1] x [0, 1] phase-space square as: 

q' = 2q- [2q] (10) 

p' = ^{p+[2q]) (11) 

where the square brackets symbolize the integer part of a real number. This evolution has a very simple geometrical 
interpretation, as a "stretching" step followed by "cutting" step, as a baker rolling dough. The map is uniformly 
hyperbolic, with a single Lyapunov exponent, with a value of In 2. Moreover, at every point the stable and unstable 
manifolds are parallel to the coordinate axes. The baker's map has a remarkably simple symbolic dynamics, and can 



2 



be mapped into an unrestricted Bernoulli shift on two symbols [ [l3[ . If the phase space coordinates p and q are written 
in binary notation 



q = O.eoeie2e3 



p = 0.e_ie_2e-3e-4 ■ 
. e_4e-3e-2e-i • eoeie2e3 ■ 



(12) 



The map's action on these symbols is to move the most significant bit of q to p, shifting to the right the decimal point 
in q: 



gr' = 0.616263 •■ • _p' = 0.6o6_i6_2e-3 ■ 

{p', q') = ... 6_4e-3e-2e-ieo • eie2e3 ■ 



(13) 



Thus each doubly infinite sequence of binary digits represents a unique trajectory. The phase space points on this 
trajectory are obtained by placing the dot somewhere (the present) and reading off the coordinate and momentum to 
the right and left of it. 

The quantization procedure for the map is not unique but follows closely semiclassical prescriptions. As originally 
formulated it used periodic boundary conditions {xq — Xp — 0) but was later modified to antiperiodic conditions 
iXq = Xp 5) The full semiclassical theory has been developed in jl^. The resulting unitary matrix for one 

step of the map is 



baker 



2N 



Gn 
Gn 



(14) 



A different approach, leading to the same result, was proposed in ||T^ (see also jlsj) interpreting the map as the 
quantization of a Bernouilli shift showing that it could be implemented by elementary gates and thus would be an 
interesting candidate to be run as an algorithm in a quantum computer. This approach is strongly based on the 
symbolic dynamics of the map and runs as follows: as the position and momentum basis are related by a Fourier 
transform, we can implement the Bernoulli shift - which shifts the most signi ficant bit of the coordinate to the most 
significant bit of momentum - using the quantum circuit shown in figure II A . 



FIG. 1. Quantum circuit which perform a Bernoulli shift, i.e. an iteration of the quantum baker's map 



The action of this circuit can be understood as follows: the qubits on the input codify in the usual way the 
eigenvalues of the position operator with the most significant bit at the bottom; after applying the split Fourier 
transform, the most significant qubit now represents the most significant bit of the eigenvalue of the momentum 
operator. As a final step, an inverse Fourier transform allows us to look at the final state in the position basis. This 
is exactly the circuit representation of the matrix Btaker ■ As is well known p9[ the Fourier matrices can be further 
decomposed into elementary gates leading to a circuit representation in terms of operations on qubits. 



B. The Harper's map 

The baker's map does not capture the full complexity of chaotic motion in Hamiltonian systems. The fact that it 
has uniform hyperbolicity and very simple manifolds are not very generic properties as the most common situation is 
that of a complex mixture of elliptic islands interspersed by chaotic regions with locally defined Lyapounov exponents. 
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To address this more general situation a different family of maps can be devised whose characteristic is to alternately 
produce "kicks" of potential or kinetic energy. The combined action of these kicks is equivalent to that of a periodic 
time dependent Hamiltonian and the resulting motion is area preserving and can mimic the full complexity of a generic 
system. A wide variety of maps on the sphere (kicked tops), on the cylinder (kicked rotors), or on the torus (Harper's 
map) have been extensively studied in the quantum chaos literature [polpl] . 

Here we have studied a specific kicked map, known as the Harper's map because of its relationship to the Harper's 
hamiltonian in solid state physics. This map acts on the unit square with periodic boundary conditions in phase-space 
and is defined by the transformation |2l|] 



q' = q — J sin(27rp) (mod 1) 
p' = p + J sin{2TTq') (modi) 



(15) 
(16) 



The behaviour of the map is determined by the real parameter 7. If 7 ^ 1 the map approaches an infinitesimal 
transformation with approximately conserved energy and its evolution is regular. If, on the other hand, 7 « 1 the 
map becomes fully chaotic. In between, the motion presents the complex mixture of regular and chaotic motion 
characteristic of most realistic systems. 

The reason that kicked maps are so popular is that they are very easy to quantize. In fact the potential and the 
kinetic kicks are respectively diagonal in the coordinate and momentum basis and therefore the full evolution consists 
of these diagonal kicks interspersed by the Fourier transformation between these two bases. One full step of the map 
is then the evolution operator U, 



U = UqG\UpGN 



For the Harper's map the two diagonal operators Uq and Up are 



Uq\qn) = e 



— -i7^'cos(2^(„+Xp)) 



qn) 



(17) 



(18) 
(19) 



and Gat is given by (g). 

Kicked maps also lead very naturally to a circuit interpretation, like in the case of the baker's map. The diagonal 
interactions are controlled phases that act among all the qubits while the Fourier transformations can again be 
decomposed into elementary gates. Some of these maps (the cat maps, for example) can be efficiently decomposed in 
terms of elementary o pera tios and have recently been studied as possible candidates to be simulated in a quantum 



computer 



Figure 



II B shows the structure of a quantum circuit which implements a generic kicked map. 



FIG. 2. Quantum circuit which implements a quantum kicked map. 



III. DISSIPATIVE MAPS 



Our aim is to study the impact of the process of decoherence induced by the interaction between our system and 
an external environment. The system would otherwise evolve according to one of the unitary operators described in 
the previous section. In general, modeling the coupling to the outside world may be complicated. Here, we will not 
use a microscopic model of this interaction but will describe the effect of the environment in a phenomenological way 
by defining a dissipative map for the evolution of the density matrix of the system. As decoherence generally induces 
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a loss of purity, the first important point to notice is that the state of the system should be defined in terms of a 
density matrix p. In the absence of any coupling to the environment, the density matrix p evolves unitarily according 
to the map 

p' = UpU'^. (20) 

The coupling to the environment would induce nonunitary evolution. To correspond to an allowed temporal evolution 
(that should come from a unitary map for the whole Universe), the nonunitary map for the density matrix has to 
satisfy some constraints. Assuming that: i) it is linear and preserves hermiticity, ii) it is trace preserving and iii) it 
preserves the complete positivity of the density matrix p9| ; the map is strongly constrained. Moreover, if one imposes 
also Markovian behavior (neglecting all memory effects) one can show that the most general map for the density 
matrix of the system should be of the form: 

p' = Y,EupEl (21) 

k 

Equation ( ^l| ) is known as the operator sum representation (or Kraus representation) of the superoperator which 
maps p onto p' — %{p) (see for a review of this representation and a derivation of the main formulae). The 

trace preserving nature of the superoperator defines a constraint for the operators Ek'. 



E 



ElEk = /, (22) 



where / is the identity operator. There is no other constraint on these operators, although it is worth noting that this 
constraint boils down to plain unitarity if there is only one operator. There is a nice physical interpretation for the 
operators: One can think equation (|2l| ) as corresponding to a process where the state p is converted randomly 

into the state fe;^^^^) with a probability tr (EkpEl\ (in this sense, these operators are quantmn jmnp operators). 

It is important to notice that the set of the Ek corresponding to a given non-unitary evolution is not unique: For 
example, if Uki are the elements of a unitary matrix, we can define the operators E'^ as 

E'l^Y^^kiEk, (23) 
k 

and show that the set ^E'^^ generates the same evolution as the set |-^;|- 

We will use the operator sum representation to define a specific model to introduce decoherence in our system. 
We want our model to correspond, in the continuum limit (where TV — > oo) to a diffusive environment having similar 
effects than the ones present in the well studied Brownian motion model |^^. For this, we will assume that the 
temporal evolution is divided in two steps: a unitary step, where the density matrix evolves unitarily as in (^); and 
a dissipative step, where the density matrix evolves by a map whose operator sum representation is of the form: 

p' = {l- a)p + ^{D{Aq, Ap)pD\Aq, Ap) + D^Aq, Ap)pb{Aq, Ap)). (24) 

where a is a real number between and 1 measuring the strength of the coupling to the environment and D{Aq, Ap) 
is the displacement operator defined in (^, which from now on will be denoted simply as D. Notice the constraint 
( p^ ) is automatically satisfied because D is unitary. In fact, in the above equation the three terms appearing in the 
operator sum representation are: 

Eo = V(l (25) 
El = (26) 

E2 = (27) 

which are normalized in such a way that (^) holds. 

The above operator sum representation provides an intuitive interpretation for the evolution: the density matrix is 
first evolved with the unitary quantum map (the unitary step). Then, three things can happen: i) with probability 
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I ~ a the state does not change; ii) with probabiHty ^ the state is displaced in one direction in phase space; iii) with 
probabihty ^ the state is displaced in the opposite direction in phase space. As the probabilities of both displacements 
are equal, a localized state does not drift in phase space as a consequence of this evolution. The net effect is to smear 
the state in phase space in the direction of the displacement operator. Thus, this nonunitary map is a discrete model 
for a diffusive process. The diffusive superoperator can be made a more efficient by using not just a single displacement 
operator but rather a sum of many terms. Thus, we will consider a more general model where 

M 

p' = (1 - a)p + ^ ^ {D^pdJ + dJpD^) (28) 

n=l 

where _D„ = D{Aqn, Ap„) for some displacements Ag„ and Ap„. In particular, it is simplest to consider all displace- 
ments along the same direction (i.e. Aqn = nAq and Ap„ = nAp). In such case, we can get the following formula for 
the evolution of the matrix elements of p (in the basis of eigenstates of 

From this equation we see that the effect on the density matrix in this representation is clear: the diagonal elements 
are unaffected by the non-unitary evolution while the non-diagonal elements are suppressed by a factor that rapidly 
decays with the distance to the diagonal. In the continuous case this suppression is gaussian but here, due to the 
discreteness of the Hilbert space the suppression has the shape of a diffraction-like kernel. Using the above expression 
it can be shown that if M — N, the net effect of one step is to completely wipe out all the non-diagonal elements, 
leaving a diagonal density matrix. 

The diffusive character of this non-unitary map is even clearer if we represent the quantum state in phase space 
using the Wigner functions. As we are dealing with a finite dimensional Hilbert space we should use a discrete version 
of the ordinary Wigner function which is well adapted to the finite phase space structure ||2^ . The definition and 
main properties are given in the Appendix and more in detail in Ref . |27| . In this representation the Wigner function 
of the density operator is a real array defined on a 2N x 2N lattice in phase space which shares many of the well 
known properties of the continuous Wigner function. The action of the diffusive map given in equation ( p8| ) in the 
Wigner representations is 

M 

W\q,p) (1 - a)W{q,p) + ^ ^[I^(g + 2nAq,p + 2nAp) + W{q - 2nAq,p- 2nAp)] (30) 

n=l 

This expression shows clearly that the diffusive step smears the Wigner function in directions specified by Ap, Aq (the 
factor of 2 present in the above equation is due to the fact that the phase space is an array of size 2A^ x 2A^ rather 
than N x N). More general diffusion models involving mixtures in different directions in phase space, or even over 
whole areas where diffusion occurs are similarly represented. 



IV. RESULTS 



In this section we present and discuss our results. They concern two separate but connected issues: First, we 
will examin the behavior of entropy as a function of time focusing both on understanding the role of the different 
mechanisms that make entropy to grow and determining the dependence of the rate at which it grows on the parameters 
of our model (like the strength of the coupling, etc). Second, we analize the issue of the correspondence between 
quantum and classical dynamics. The basic tool we used for our studies is a highly efficient code to evolve the density 
matrix which makes good use of fast Fourier transform routines. With this, running on a PC, one can easily compute 
the density matrix for Hilbert spaces with dimension of about N — 2000 (larger machines would be required to go 
above that limit but there is no good reason to do that, see below). The source code is available from the authors. 



A. Rate of entropy production 



For a variety of initial states we studied the evolution of the density matrix under the dissipative map and computed 
the linear entropy S = — \og{Tr{p^)). We did that both for the baker's and the Harper's map. The diffusion mechanism 
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is as in ( p8|) and in the plots ap labels diffusion along the momentum axis {D — V ) while aq labels diffusion along 
the coordinate {D = U) Our results clearly establish that in both cases, the rate of entropy production becomes 
independent of the parameter a provided its value is above a certain threshold (that is close to 0.4). This is seen 
in Figures ^ where the plots show the entropy as a function of time for two types of diffusive environments (where 
diffusion is along either position or momentum). It stands out from this graph that all of the curves with a greater 
than 0.4 have a linear regime with the same slope. This slope turns out to be equal to In 2, the Lyapunov exponent 
for the baker's map. 
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FIG. 3. Linear entropy growth for the baker's map with N=1594 and different values for a. For the top graph diffusion is 
along momentum and for the bottom graph along position. 

Similar results are shown in Figure ^, for which the unitary propagation was provided by the Harper map. Even 
though in this case there is not a well defined linear regime (at least not so well defined as in the case of the baker's 
map), it is clear that the rate of entropy production becomes independent of a. Thus, in this regime entropy is 
produced due to the coupling with the environment but the rate becomes independent of the coupling strength. 




7 



7 



- - a =0.1 

p 

_ _ a =0.3 
p 

- - a =0.5 

p 

a =0.8 

p 

a =1 



25 




FIG. 4. Same as Fig H but for Harper's map (7 — 0.45) 



This is one of the main resuhs of the paper, that substanciates the conjecture that was originaUy put forward in 
1^ and later tested numericaUy in various works (|J^;|l^ • 

We have also studied the dependence of the entropy growth as a function of N (the dimensionality of the Hilbert 
space) . This was done both for small values of a and in the regime where the rate of entropy production is independent 
of the coupling. These results are shown in Figures I and H In both cases we have found that - when appropriately 
rescaled - the entropy growth is quite independent on the value of N and tends to a universal curve in the limit 
N ^ 00 (as 1/N is the effective Planck constant this corresponds to the classical limit). This result suggests that the 
mechanism for entropy growth in quantum maps has a dominant classical component with corrections that vanish in 
the classical limit. 
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FIG. 5. Scaling of the linear entropy growth with dimension of the Hilbert space of the system with diffusion on the position 
axis. The strength of the diffusion is kept constant by smoothing over a constant fraction of the system's phase space. In 
relationship with Figures ^ and ^ both axis have been rescaled by log(A'^) 
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FIG. 6. Scaling of the linear entropy growth with dimension of the Hilbert space of the system with diffusion on the 
momentum axis. The strength of the diffusion is kept constant by smoothing over a constant fraction of the system's phase 
space. In relationship with Figures ^ and |^ both axis have been rescaled by log(Ai'). 

Our results show clearly that above some threshold for a, the rate of entropy production is determined by the 
dynamics of the unitary quantum map, and in great measure by its classical counterpart). Why is this the case? We 
will argue here that it is possible to attribute the entropy production to two interconnected processes whose origin 
can be better explained by using a phase space representation for the quantum state. The reason is that in this way 
we can use some of the intuition we have about the behavior of the classical system. Consider an initial state that 
is represented by a localized and smooth phase space distribution (for concreteness we base our discussion on the 
Wigner function, see below). The application of the chaotic unitary map will distort the state in a way that, at least 
for short times (times smaller than log(7V)), will be consistent with the classical evolution. In a hyperbolic region this 
will mean that the initial wavepacket will be stretched along the unstable manifold and contracted along the stable 
one. In the case of the baker's map (and only in this case) these two directions are the same all over phase space: the 
vertical direction -that we have called momentum - is stable (contracting) and the horizontal one - the coordinate - is 
unstable (expanding). The rate of expansion is equal to the Lyapunov exponent which in the baker's map (and only in 
this case) is constant all over phase space and is equal to In 2. However, this expansion cannot proceed forever. As the 
phase space is finite, after some time the state will start experiencing the other essential feature of chaotic dynamics: 
folding. As soon as this process starts the quantum and classical evolutions will start deviating significantly from 
each other due to the quantum interference between the different folds. In the Wigner function (see below) these 
interferences appear as rapid oscillations with large negative values happening in between the folds and oscillating 
in directions parallel to them. These oscillations are the footprint of the quantum interference between the different 
pieces of the delocalized state. 

The two interconnected effects generated by the diffusive map are: (a) the destruction of the interference fringes in 
the Wigner function, (b) the smoothing and widening of the regions where the Wigner function is already positive. It 
should be clear that the two effects have the same origin: diffusion, but it is worth distinguishing between them because 
they are effective for different types of states and can be analyzed separately under some special circumstances. In 
fact, process (a) affects states which have a quantum nature exhibiting important oscillations in the Wigner function. 
On the other hand, process (b) is present even for classical systems (in fact, in a classical difussive map process (b) is 
the only source of entropy). The important point is that even if we consider an initial state where only process (b) is 
important (i.e., a classical state), the dynamics of the system will create large scale quantum interference which will 
trigger process (a). 

In the baker's map the two processes can be analyzed separately. Thus, if we consider a non-unitary map producing 
diffusion only along the momentum (i.e., Aa; = in (p8|)) only process (b) will be relevant. In fact, this is the case 
because an initial phase space distribution will tend to be contracted in the momentum direction while being expanded 
in the position direction. Folding will start to occur when the state reaches the boundary of the phase space. After 
that, quantum interference will develop between the different pieces of the packet and will be characterized by fringes 
which will be aligned along the vertical direction. Therefore, the fringes will not be significantly affected by diffusion. 
On the contrary, the effect of diffusion will be to fight against the contraction of the wavepacket along the momentum 
direction. One can see that diffusion will tend to compensate contraction and the Wigner function will tend to 
acquire a critical width along momentum. When this happens, the entropy will start growing at a rate fixed by the 
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expansion rate which is the Lyapunov exponent. This is the scenario originally discussed in ^J: the entropy grows at 
the rate fixed by the Lyapunov exponent when the Wigner function reaches a critical width along the stable direction 
(contraction stops but expansion continues almost unafected by diffusion, therefore area grows exponentially which 
implies that entropy grows linearly). 

The other extreme case is to consider a non-unitary map where difussion is along the position direction (i.e., Ap = 
in (p8|). In that case, diffusion will not prevent the system from contracting along the stable direction and therefore 
will not give rise to a critical width in that direction. However, entropy will still grow but it's origin is entirely due to 
process (a). In fact, as soon as the phase space distribution starts to fold, interference fringes develop. These fringes 
are aligned along the momentum axis (since they correspond, roughly, to interference between two horizontal strips 
separated by a distance 1/2). As the strips are macroscopically separated they will be very sensitive to diffusion along 
position. Thus, if the coupling strenght a is such that the fringes are washed out in just one iteration the mechanism 
will produce a mixture of two horizontal strips every iteration of the map producing one bit of entropy. The rate 
associated to this mechanism is the rate of folding that in the case of the baker's map (and only in this case) is equal to 
the global Lyapunov exponent. Thus, the fact that both processes (a) and (b) give rise to the same rate is an accident 
of the baker's map. In general the rate will not be the Lyapounov exponent but a combination of such exponent and 
the global folding rate. In any case, the general conclusion which is still valid is that the entropy production rate for 
classically chaotic systems will tend to be dominated by dynamical properties of the system becoming independent of 
the strength of the coupling to the environment. 

It is interesting to present a simple analytic model to support the above discussion about entropy production. The 
simplest case is the one corresponding to process (b), where difussion is along the position direction. To simplify 
our model we will assume that the action of the baker's map is just to transform a momentum eigenstate into a 
superposition of two such states: 

^|M = ^(I§> + I^» (31) 

(this is a rough approximation for the baker's map that neglects diffraction effects and relative phases but that 
captures the essential folding mechanism of the baker's action in the large N limit and away from saturation ). In this 
approximation, the unitary evolution turns each momentum state into a coherent superposition of two momentum 
states separated by a macroscopic distance iV/2. According to (p9[), the momentum diffusion leaves untouched the 
diagonal - incoherent - elements for the density matrix of this state while suppressing, in the large M limit, the 
off-diagonal parts by a factor 



With these assumptions it is simple to obtain an equation for the entropy as a function of time. Thus, after t 
iterations the density matrix has 2*+"^^ matrix elements with value 2~*(1 — a)", with n ranging from 1 to t. Using 
this, we showed that the linear entropy is 

S'(t) ==tln2-ln[(l-a)2(l + (2(l-a)2)*)-l]. (33) 

The second term in this equation corrects the slope in the small a regime. The transition is predicted at a critical 
value of Q!c = 1 — 2~^/^ w 0.3. For a > the second term becomes intependent of t and the slope saturates at 

= ln2. In the weak coupling regime , for a < etc the growth is still linear, but the slope depends on the coupling 
strenght and is given - for long times - by 

S' = -21n(l-a) (34) 

In this regime the slope of the entropy is then linear with a - and independent of the Lyapounov exponent of 
the map. All the predictions of this simple model are confirmed by the numerical calculations. In Figure (0) we 
display the behavior of the rate S* as a function of the coupling strength a during the regime of linear growth. The 
agreement between the simple model and the exact numerical results is very good, describing both coupling regimes 
very accurately. To see the limitations of the model we plot in Figure (|^ the same quantities, but taken at later times, 
closer to the point where the entropy saturates due to the finiteness of the model. Here the agreement deteriorates 
as the assumptions of the model break down. 
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FIG. 7. 5*2 vs. a: Analytic prediction against numerical results. The analytic line and the numerical data points were not 
taken from the same iteration, in order to compensate for the fact that the initial state in the simulations was not a momentum 
eigenstate but a coherent state instead, which takes some iteratios to transform into something approximate to a momentum 
eigenstate. 




X X X 



Theoretical prediction 

X Numerical data 



0.4 ^ 0.6 

q 



FIG. 8. 5*2 vs. a: Anaylitic prediction against numerical results some iterations after the ones depicted in graph |^ 



B. Quantum classical correspondence. Wigner function 

We also studied the evolution of the Wigner function for a variety of initial states both for the baker's map and the 
Harper's map. In Figures 0and|ll| we compare the classical (left), unitary (center) and quantum difussive (right) 
cases for these two maps. The general features that are observed are the following: The unitary map follows the 
classical evolution for up to a time which is of the order of log (TV). After this time quantum interference develops 
between the different pieces of the Wigner function. These effects are responsible for the loss of correspondence between 
quantum and classical evolution that was discussed in H,^-^ . When the effect of the coupling to the environment 
is taken into account, it is clear that the correspondence between the quantum and the classical evolution is restored 
and holds for a much longer time. However, it is worth pointing out that the fact that our system is finite (i.e., 
it has N orthogonal states) implies that saturation is achieved and the system approaches an equilibrium state. 
This equilibration takes place for times which are a few times log(iV) as seen in the curves for the entropy (thus, 
as entropy grows linearly, it will approach the log(iV) value relatively soon). Therefore, in this type of system the 
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quantum classical correspondence is restored but has a rather short interesting dynamical regime. The nature of the 
final equilibrium state has footprints of the underlying classical dynamical system that are evident in the presence 
of decoherence. This is not so transparent by analyzing the final state for the Baker map but it is evident for the 
Harper's map: the final state with decoherence is an equilibrium state where occuping uniformily all the available 
phase space with the remarkable exceptions of the regular islands. 

By analyzing the Wigner distribution we can study the differences between the two processes contributing to the 
entropy growth that were discussed above. Thus, when difussion is along the direction of momentum it tends to 
prevent the contraction of the Wigner function but does not affect the interference fringes. On the other hand, when 
diffusion is along the position direction it does not affect the contraction but efficiently destroys the interference 
fringes. In Figure ^ we show a comparison of the Wigner function for the cases when diffusion is in position and 
momentum. It is clear that in the first case the width in momentum is significantly smaller than in the second while 
the interference fringes (that are alligned along momentum) are washed out much more efficiently. When the diffusion 
is along the stable manifold, the packet width does not shrink in momentum and approaches a critical width. Besides, 
the interference fringes are more noticeable. In the case of position diffusion there is also a critical width that in this 
case corresponds to the wavelength of the interference fringes that are no longer efficiently suppressed by diffusion 
(that damps all small wavelength fringes). 



>4lll!li!l; 



Position diffusion 



Momentum diffusion 



FIG. 9. Comparison between diffusion in the different manifolds of the baker's map. 



V. CONCLUSIONS 



The main results in this paper come from numerical simulations of open quantum maps, particularly the baker's and 
Harper's map coupled to a diffusive environment. Our simulations provide strong numerical evidence supporting the 
conjecture put forward in when a quantum system whose classical analog is chaotic is subject to decoherence, there 
is a regime in which the rate of entropy production is determined by the system's dynamical properties (Lyapunov 
exponents, folding rates, etc) and becomes independent of the strength of the coupling to the environment. Another 
interesting aspect of our results is the existence of a clear scaling for the entropy as a function of the dimensionality 
(or of the effective Planck constant). The scaled curves for the entropy become independent oi h = jj) with quantum 
corrections that decay as N grows. We also found cualitative evidence, looking at the phase space distributions, that 
the introduction of diffusion extends the time of coincidence between the classical and quantum evolutions. Our 
results also clarify the origin of entropy production allowing us to develop new intuitive explanations to understand 
the origin of the entropy growth (see pl|] for another point of view based on the use of consistent histories). The 
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two processes that are responsible for such growth (washing out of interference fringes and smoothing of the positive 
pieces of the Wigner function) can be studied separately only in the case of the baker's map where stable and unstable 
manifolds have a global nature in phase space. The differences between the behavior of the entropy as a function of 
time for both non-unitary models (diffusion along position or along momentum) arc illuminating: in both cases there 
is a regime in which the entropy grows linearly with a slope that is equal to the Lyapunov exponent (which is also 
the folding rate). The initial behavior is different in the two cases: when diffusion is along the position direction the 
entropy stays small initially until the wavepacket reaches the boundary of phase space and folding becomes effective. 
Then every iteration of the map creates new interference fringes that are rapidly destroyed by decoherence. Thus, the 
dominant process for entropy production is (a). When diffusion is along the momentum direction, entropy growth 
starts at early times. In this case the growth of entropy is produced by process (b) and is associated with the fact 
that diffusion tends to balance the effect of the contraction of the wavepacket along the stable direction giving rise 
to a critical width. In any case, our study shows that entropy growth at a rate fixed by the dynamics of the system 
rather than the coupling strength to the environment is, as first conjectured in ||^ a characteristic of quantum chaos. 
The differences between this behavior and the one that characterizes integrable systems (or quantum states localized 
in stability islands like the ones present in Harper's map) are quite dramatic (see | |3^ for related results). 

This work was partially supported by grants from Ubacyt, Anpcyt, Conicet and Fundacion Antorchas. JPP 
gratefully acknowledges the hospitality of ITP Santa Barbara where the final version of this manuscript was completed. 



APPENDIX A: DISCRETE WIGNER FUNCTION 

We use a discrete version of the Wigner function, which we describe in more detail in ||27|] . Here we give the basic 
definitions to make the paper self-contained. The Wigner function is usually defined, for a continuous system, as 

^^'^'^^^ 2^/_r^''"2'^''^+2^''*''^^- ^^^^ 

This function is the expectation value of the so-called phase space point operators: 

W{q,p)^tiipA{q,p)). (A2) 

The operators A{q,p) have a number of interesting properties from which the properties of the Wigner function can 
be derived (in fact, they are hermitian operators, they form a complete basis, etc; see |p7| ). It turns out that in order 
to generalize all these properties to the discrete case it is necessary to use a phase space lattice of 2A^ x 2N points 
and to define the phase space operators (in terms of the shift operators used in the body of the paper and defined in 
(i) as 

A{q,p)^^WRV~Pe^pi7rqp/N. (A3) 

Using these discrete phase space point operators the Wigner function has all the desired properties (it is real, it can be 
used to compute expectation values of operators and gives rise to all the correct marginal distributions when summed 
over any set of lines in the phase space grid). To explicitly compute this function one can use a variety of formulae 
among which we found convenient the following expression: 

W{qJ,p^) = ^e-^^^^-^^^"'^') X 

n fc 

On the other hand, it takes a atraightforward calculation to show that the A{q^p) operators transform in the following 
way when subjected to translations (both in their continuous and discrete versions): 

D{Aq,Ap)M<l,P)Dl^g,Ap) = Mq + 2Ag,p + 2Ap) (A5) 



We can use this property and the definition of the Wigner function given in equation A2 to derive the equation ^ we 
use in the paper. 
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FIG. 10. Phase space distribution at several steps in the evolution of the Baker map. The quantum versions have been 
computed with = 62. 
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FIG. 11. Phase space distribution at severals step in the evolution of the Harper's map. The quantum version has diflusion 
parallel to the position axis and N — 62, Uq — 0.5, M — 2. 
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